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BACKGROUND 


Early  in  1977,  Corputer  Sciences  Corporation  (CSC)  conducted  a  study 
for  the  Life  Support  Engineering  Division  (LSED)  (603) ,  Naveil  Air 
Developirient  Center  (NMX^)  under  Contract  N62269-75-0001,  Task  Order 
No.  47.  The  purpose  of  the  stut^  was  to  examine  the  feasibility  of  vising 
a  microprocessor  based  system  to  control  the  functions  of  a  vertical 
seeking  ejection  system.  In  November,  1977,  CSC  issued  a  report  containing 
the  results  of  that  study,  titled  "Microocrcputer  Controlled  Ejection 
Seat  Feasibility  Study".  CSC  found  the  caicept  to  be  feasible  and 
reccmmended  further  analysis.  To  this  end,  NAIX!  initiated  the  purchase 
of  the  IMSAI  PCS  80/30  microprocessor  based  software  develqpnent  system. 

From  January  ,  1978  to  July,  1978,  CSC,  under  Ccaitract  N62269-78-C-0191 
Task  Order  No.  7,  conducted  a  second  effort  for  NADC  and  in  July,  issued 
a  report  entitled  "Esc^ie  System  Trajectory  Simulation  and  Microprocessor 
Control  ^stem".  This  report  contains  the  result  of  that  effort,  the 
purpose  of  vbich  was  to  define  the  functional  requirements  of  an  ejection 
seat  system,  examine  the  reliability  and  maintainability  aspects  of  the 
INTEL  8085  microprocessor  for  military  use  and  outline  a  software 
demonstration  program  to  perform  the  vertical  seeking  maneuver  and  to 
control  the  timing  and  sequencing  of  specific  events  during  ejection. 

The  data  used  to  demonstrate  the  program  was  provided  by  the  Naval 
Vfeapon  Center,  China  Lake,  CA. 


( 


In  Sepbeniber  1978,  CSC  corpleted  the  developnent,  testing,  and  demonstra¬ 
tion  of  the  program.  In  October,  1978,  CSC  issued  an  interim  technical 
report  entitled  "Supplement  to  Esca^  System  Trajectory  Simulation  and 
Microprocessor  Control  Systan",  which  documents  the  program  that  was 
developed  and  tested  to  demonstrate  the  capabiliiy  of  the  INTEL  8085 
microprocessor  to  control  the  functions  of  a  vertical  seeking  ejection 
seat. 

In  June  1979,  CSC  was  asked  to  implement  the  simulation  of  the  vertical 
seeking  ejection  seat  in  tiie  existing  Aircrew  Escape  Systan  Conputer 
Simulation  Program  -  ICARUS.  In  October  1979,  CSC  ccmpleted  the 
implen^tation  an!  preliminary  testing  of  the  revised  version  of  the 
ICARUS  Program. 
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1.0  DJTHXmcriQN 

Three  major  functions  were  added  to  the  ICARUS  Program  to  give  it  the 
capability  of  simulating  the  vertical  seeking  ejection  seat;  a  system 
viiich  continuously  measures  the  angular  rotation  of  the  seat/occupant; 
a  procedure  to  update  the  seat/occvpant  position  through  a  direction 
cosine  update  algorithm;  a  oormiand  decision  function  vdiich  determines 
the  rocket  roll  and/or  pitch  commands  ty  using  the  values  in  the 
direction  cosine  matrix  and  the  angular  rates  of  the  seat/occvpant 
system.  Ihe  displacenent  of  the  rocket  is  required  to  get  the  seat/ 
occi:5)ant  system  into  an  upright  position  within  the  shortest  time 
possible. 

The  material  presented  herein  assxmes  that  the  reader  is  familiar  with 
the  input  file  and  basic  structure  of  the  ICAHJS  Pix^gram. 
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2.0  VERTICa^L  SEEKING  MANEUVER 

The  Vertical  Seeking  Algorithm 
5  basic  steps  cotprising  the  vertical  seeking  algorithm. 
Initialize  Seat/Occupant  Orientation. 

Integrate  Rotational  Rate  Data. 

Update  Direction  Cosines. 

Generate  Rocket  Giraballing  Conmands. 

Gimbal  Rocket 

2.2  Initialize  the  Orientation  of  the  Seat/Occupant  System 

Two  coordinate  systems  (shewn  in  Figure  1)  are  involved  in  the  computation 
of  the  vertical  seeking  ejection  system,  the  Seat/Man  Coordinate  System 
(SM:::S)  and  the  Vertical  Seeking  Coordinate  System  (VSCS) .  There  are  three 
elements  of  the  direction  cosine  matrix,  D3(l) ,  D3(2) ,  and  D3(3) ,  vdiich 
are  defined  at  ejection  initiation. 

The  direction  cosine  matrix  (DVC)  which  transfers  the  Z-Axis  of  VSCS  to 

the  SbCS  is  given  by 

Cos(C7)  0  Sin(C7) 

0  10 
-Sin(C7)  0  Cos(C7) 

The  direction  cosine  matrix  ASEM  viiich  transfers  a  vector  from  the  VSCS 

to  the  Earth  Fixed  Coordinate  System  (EPCS)  is  given  by 

[asem]  =  [dvc]  X  [ase] 


DVC  = 


2.1 

There  are 

1. 

2. 

3. 

4. 

5. 
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SMCS  -  SEAT/M?^  COORDINATE  SYSTEM. 

VSCS  -  VER3TCAL  SEEKING  aXlRDINATE  SYSTEM. 

~  ANGTiE  (\^1RT/  ^'\/SCS^ 
C^2  =  ^^NGLE  (vert/  Yyg^) 
I  ^33  ”  ANGEiE  (VERTf  ^VSCS^ 


Figure  1  -  VSCS,  SMCS  Definition 
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ASE  is  a  direction  cosine  matrix  viiich  transfers  a  vector  fran  SMCS  to 
the  EPCS.  The  initial  values  of  the  three  direction  cosine  matrix 
elements,  D3(l),  D3(2),  and  D3(3)  are  calculated  as  follows: 

3  2 

D3(I)  =  -ASEM  (1,2)/ (  >  ASM  (I,J)  )  ^ 

J=1 

Where  I  =  1,3. 

2.3  Integrate  Rotational  Rate  Data 

The  rotational  rates  cire  integrated  with  respect  to  time  in  order  to 
keep  track  of  the  rotation  of  the  seat/occupant  syston  about  each  of 
the  axes.  The  integration  scheme  is  based  on  the  following  formula 
for  the  trapezoidal  rule:  •  ' 

A  =  %  AhCYj^  +  2Y2  +  2Y2  +  •  •  •  •  +  Yjj) 

where  Y^  corresponds  to  rotational  rates  for  a  particular  tine.  Since 
data  is  sampled  approximately  every  1.25  ms, Ah  =  1.25. 

The  piogram  does  not  update  the  rotation  position  about  each  of  the  seat/ 
occupant  axes  at  every  time  increment  but  only  vdien  the  rotation  has 
exceeded  a  predefined  angle.  A  is  defined  as  the  maximum  angle  that 
could  be  generated  in  one  sampling  period  of  rotating  at  the  maximum 
rate  that  can  be  sensed  by  the  gyroscopes.  In  this  case,  the  maximum 
rate  that  can  be  sensed  by  the  gyroscopes.  In  this  case,  the  maximum 
rate  is  8.72  rad/sec  and  the  sampling  period  is  1.25  ms;  therefore 
A  =  .0109  radians.  If  hovrever,  the  seat/ocCupant  system  is  rotating  at 
the  maximum  rate  in  a  negative  direction  (-8.72  rad/sec) ,  then  A  =  -.0109 
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radian.  When  the  seat/occxpant  system  has  rotated  at  least  .0109  rad. 
in  either  direction  around  any  of  the  3  axes,  the  appropriate  direction 
cosines  are  updated  and  the  corresponding  test  value  is  subtracted  from 
the  sum;  the  integration  then  continues  as  before. 

2.4  Update  Direction  Cosines 

Implementation  of  this  part  of  the  vertical  seeking  algorithm  is  based 
on  the  CrowdeT’-Hession  Direction  Cosine  l^xJating  Algorithm  as  follows; 

^31, a  =  C31,k 

^32, a  =  ^32,k  “  ^31, k 

^33, a  =  ^33, k  +  ^31^^ 

=  C3i^a  +  A02  C32,a 

C32,b  =  C32,a 
C33,b  =  C33,a  -  A0x  C32,a 
C31,k+1  =  C3i^b  -  AGy  C33^j3 

C32,k+1  =  C32,b  +  C33,b 

C33,k+1  =  C33,b 

Where; 

^31,  ^32,  ^33  three  direction  cosines  D3(l) ,  D3(2) ,  D3(3) 

K,  K  +  1  designate  points  in  time; 

a,  b  are  values  for  intermediate  calculations  used  in  arriving  at  the 
next  point  in  time  and 

^®X,  angular  displacenents  of  the  gyros. 

The  displacements  are  represented  by  a  quantum  pulse  weight  A0  which 
has  the  general  form  of  2“^,  vhere  n  is  determined  by  the  saitpling  time 
and  the  maximum  value  of  the  angular  rate  is  given  as  7  (2“^  =  .0078  rad) . 
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Every  2.5  itis,  this  i^date  algorithm  is  processed  twice.  In  processing 
the  algorithm,  A0  ^  and  ^  are  determined  by  the  integration  of  the 

rotational  rates  around  the  X,  Y,  and  Z  axes  respectively.  If  the  seat/ 
occupant  systan  has  rotated  less  than  .0109  radians  in  either  the  positive 
or  the  negative  direction,  the  corresponding  A0j  =0.0.  If  the  seat/ 
occ\:5>ant  system  has  rotated  at  least  +  0,0109  radians,  the  corresponding 
A0j  =  +  2~^  -  +  0.0078  radian. 

2.5  Generate  Rocket  Gdmballinq  Conmands 

TXfTO  oomtBnds,  a  roll  conmand  PC  and  a  pitch  conmand  QC,  are  output  to 
ginibcil  the  rocket  and  are  governed  by  a  set  of  rules  as  follows; 

1.  If  the  Z-Axis  of  the  VSCS  is  parallel  to  the  normal  to  the 
earth  surface,  a  maximum  roll  ccmmand  of  2.0  is  generated; 

2.  If  the  seat/occupant  is  in  an  inverted  position,  (D3"(3) ,  <  0.0) , 
a  large  ccmmand  (between  1  and  2)  is  generated  in  order  to 
rotate  the  seat/occupant  back  into  the  upright  position  i 
quickly; 

3.  If  the  seat/occupant  is  in  an  U5)right  attitude,  a  moderate 
pitch/ioll  ccmmand  is  generated  according  to  the  values  of 
D3 (1)  and  D3 (2) ; 

4.  The  extreme  values  for  either  PC  or  QC  are  limited  to  +2.0 
and  -2.0. 

The  rocket  gimballing  coiinand  control  law  is  given  in  Figure  2. 
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|D3(1)  1  |D3(2)  I 

PC  =  -D3 (2) 

QC  =  (2-|D3{1) 1) 

*  SIGN  (D3(1)) 


Figure  2  -  Rocket  Gimballing  Cotmand  Control  .  Law 
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Giitibal  Rocket 


The  seat-rocket  motor  arrangement  as  inplemented  at  MC  is  shown  in 
Figure  3.  It  can  be  seen  from  the  figure  that  the  rocket  motor  can 
be  girtiballed  +  16°  in  two  directions. 

The  roll  ccmraand,  PC,  and  the  pitch  ccramand,  QC,  are  output  to  a 
Signal  converter,  vMdi  in  turn  causes  the  rocket  to  move  along 
the  X-Axis  and/or  Y-Axis  of  VSCS.  The  displacement  of  the  angle  is 
measured  not  from  the  previous  position  of  the  rocket  but  from  the 
neutral  position. 

# 

There  are  two  sets  of  equations  involved  in  simulating  this  signal 
conversion; 

The  roll  thrust  vector  control  angle,  X14,  is  calculated  as  follows; 

~  ^^76  *  ”  ^^77  *  ^4^ 

%6  =  T  ■  *  ’'l6> 

fr  =  5  *  +  (  20  ‘  Xj^g)) 

X35  =  (31.45  *  -  (31.45  ♦  X35) 

^36  '  ’'35  ■  *  =^4>  ■  *  ^36* 

^14  “  ^36 

is  obtained  tlirough  integration  over  X^^. 

X4  is  the  bo<^  roll  rate  and  PC  is  the  roll  ccramand. 

The  pitch  thrust  vector  control  angle,  calculated  as  follows: 


-10- 


SEAT  iVIOTOR  /.  SRAWGEIVI EIMT 


©r  =  (C74  *  QC)  -  (C75  *  X5) 

45  =  ®t-/ioo*V 

%  =  5  *  +  (20  *  X^g) ) 

X40  =  (31.45  *  Qj,)  -  (31.45  *  X4Q) 

^41  =  ^40  "  (50992  *  X3^3)  -  (48.26  *  X4j^) 

^13  =  ''41 

The  pitch  angle  is  obtained  through  integration  over  X13.  X5 
is  the  boc^  pitch  rate  and  QC  is  the  pitch  coimand. 
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3.0 


PROGRAM  DESCRIPTION 


The  following  disciassion  describes  the  modifications  that  were  made  to 
the  ICARUS  program  in  order  to  inclvjde  a  vertical  seeking  maneuver 
simulation  capability. 

3.1  Program  ICARUS 

common  blocks  were  added  to  program  ICMIUS  in  order  to  provide 
the  inter-subTOutine  data  camiunication  related  to  the  vertical 
seeking  maneuver  simulation.  Cannon  block  MPESIN  contains  the  irput 
data  and  the  related  constants  and  cctnmon  block  MPES  contains 
intermediate  variables  used  in  corputations. 

COMMDN/MPESIN/PSM,  QSM,  RSM,  IMPES,  IMRKT,  BOT,  IMREl,  IMRE2,  C5, 

C6,  Cl,  C74,  C75,  CIS,  C77,  DVC(3,3). 

CCMOJ/MPES/D3  (3) ,  DALPH(3),'  DTH(3,2),  ANGR(3) ,  PC,  QC,  X13,  X14,  Xl5, 

X16,  X35,  X36,  X40,  X41,  XD13,  XD14,  XD15,  XD16,  XD35,  DX36, 
XD40,  XD41,  PTIME,  lUPDATE,  PPH,  PTHAT. 

3.2  Subroutine  EXBCUT 

Flags  IMPES,  IMRKT  and  IMCAT  are  set  according  to  appropriate 
irput  values. 
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IMPES  -  Vertical  seeking  seeking  maneuver  flag,  set  by  INT(IO) 
=  0  i  OFF 

=  1  ;  CN,  with  non-filter  version 
=  2  ;  QJ,  with  filter  version 

IMRKT  -  kxket  Ignition  time,  set  by  INT(ll) . 

=  1  ;  time  equal  to  TRI 
=  2  ;  calculated  rocket  ignition  time 

IMCAT  -  Catapult  flag,  set  by  INT  (12) . 

=  ;  no  catapult 

=  1  ;  catapult 

IMREl  -  Vfork  array  index  before  drogue  chute  projection. 

TMRF.2  -  Work  array  index  after  drogue  chute  projection. 

3.3  Subroutine  INIT 

A  set  of  constants  was  defined  through  data  stetenents. 

C5  =  .2792 
C6  =  .2792 
C74  =  .54 
C75  =  .09 
C76  =  .54 
C77  =  .09 
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Angle  C7,  which  is  the  angle  between  the  Z-Axis  of  the  SMCS  and  the 
Z-Axis  of  the  Vertical  Seeking  Coordinate  System  (VSCS) ,  is  defined  in 
input.  Matrix  DVC  transfers  the  data  from  the  VSCS  to  the  EMSS. 
DVC(1,1)  =  CX)S  (C7) 

DVC(1,2)  =  0.0 
DVC  (1,3)  =  SIN  (C7) 

DVC(2,1)  =  0.0 
DVC(2,2)  =  1.0 
DVC(2,3)  =0.0 
DVC (3,1)  =  -SIN(C7) 

DVC(3,2)  =  0.0 
DVC(3,3)  =  COS(C7) 

Matrix  ASE,  vhich  transfers  the  data  frcm  SMCS  to  EPCS,  along  with 
Matrix  DVC,  was  used  in  the  cotputation  of  the  initial  values  D3  (1) , 
and  D3(2)  and  D3(3)  of  the  direction  cosine  matrix. 

ASEM  (1,1)  =  (DVC(1,1)  *  ASE  (1,1))  +  (DVC  (1,2)  *  ASE  (2,1))  + 
(DVC  (1,3)  *  ASE  (3,1)) 

ASEM  (1,2)  =  (DVC(1,1)  *  ASE  (1,2))  +  (DVC  (1,2)  *  ASE  (2,2))  + 
(DVC  (1,3)  *  ASE  (3,2)) 

ASEM  (1,3)  =  (DVC  (1,1)  *  ASE  (1,3))+  (DVC  (1,2)  *  ASE  (2,3))  + 
(DVC  (1,3)  *  ASE  (3,3)) 

ASEM  (2,1)  =  (DVC  (2,1)  *  ASE  (1,1))  +  (DVC  (2,2)  *  ASE  (2,1))  + 
(DVC  (2,3)  *  ASE  (3,1)) 

ASEM  (2,2)  =  (DVC  (2,1)  *  ASE  (1,2))  +  (DVC  (2,2)  *  ASE  (2,2))  + 
^  (DVC  (2,3)  *  ASE  (3,2)) 
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ASEM  (2,3)  =(DVC  (2,1)  *  ASE  (1,3))  +  (DVC  (2,2)  *  ASE  (2,3))  + 

(DVC  (2,3)  *  ASE  (3,3)) 

ASEM  (3,1)  =(DVC  (3,1)  *  ASE  (1,1))  +  (DVC  (3,2)  *  ASE  (2,1))  + 

{me  (3,3)  *  ASE  (3,1)) 

ASEM  (3,2)  =(DVC  (3,1)  *  ASE  (1,2))  +  (DVC  (3,2)  *  ASE  (2,2))+ 

(DVC  (3,3)  *  ASE  (3,2)) 

ASEM  (3,3)  =(DVC  (3,1)  *  ASE  (1,3))  +  (DVC  (3,2)  *  ASE  (2,3))  + 

(DVC  (3,3)  *  ASE  (3,3))  _ 

D3  (1)  =  (-ASEM (1,2)) //(^EM (1 , 1) ) \  (ASEM (1,2))^+  (ASEM (1,3)')^ 
D3(2)=(-ASEM(2,2)/\/ (ASEM (2,1))^+  (ASEM ( 2,  2) )  (ASEM (2,3))^ 
D3(3)=(-ASEM(3,2))// (ASEM(3,1))2+(aseM(3,2)  )^+(ASEM(3,  3)  )2 

3.4  Subroutine  INPUT 

Angle  C7,  vMch  is  the  angle  between  the  Z-Axis  of  the  SMCS  and  the  Z-Axis 
of  the  vertical  sedcing  coordinate  ^steoa  (VSCS) ,  is  defined  by  input 
INT(13) .  Rxdcet  thrust  tables  TER  and  TLR  were  converted  fron  the  standard 
bum  time  into  the  TRBO  burn  time  for  each  individual  test  case.  The 
rocket  thrust  was  adjusted  according  to  the  difference  in  rocket  bum 
time  in  order  to  maintain  the  same  total  mcmentum  over  the  viiole  time 
period. 


The  total  rooraentum  fron  the  standard  rocket  thrust  table  is  calculated  as 


DID 


NPTS-1 


1=1 


(YrAB{I,2)  +  YTRB((I  +  1),  2))  *  (YTAB((I  +  !),!)- 
me  (I,l))/2.0 
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Where  YTAB  (1,2)  =  Rocket  thrust  at  point  I 

YTAB  (1,1)  =  Rocket  bum  time  at  point  I 
NPTS  =  Total  number  of  points 

The  time  ratio  between  the  TRBO  and  the  standard  rocket  bum  out  time 
is  R. 

R  =  TRBO/(YEAB  (NPTS,  1)  -  'mB  (1,1)) 

The  converted  rocket  thrust  and  corresponding  time  at  point  I  is  calculated 
as; 

TNEW  (1,2)  =  YTAB  (I,2)/R 

TNEW  (1,1)  =  TNEW  ( Cl-1) ,  D  +  (1,1)  -  YEAB  ((I-l),  1))  *  R 

The  total  momentum  frcm  the  converted  rocket  thrust  table  is  calculated 
as: 

NEW  =  '^^  (TNEW  ((I-l),2)  +  TNEW  (1,2))  *  (TNEW  (1,1)  -  TNEW  ((I-l),  l))/2.0 
1=2 


If  the  total  itcxnentum  NEW  is  beyond  (1  +  .02)  of  the  total  mcroentum  OID, 
the  array  TNEW  and  NEW  are  recalculated  vdth 
R.  =  R  *  NEW/OID 

Both  tie  standard  thrust  table  aid  the  converted  thrust  table  are  printed 
for  user’s  reference. 
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3.5 


Subroutine  ROCKFM 


a)  5|ie  angular  increment  for  every  time  increment  of  0.00125  second 
is  oonputed  as; 

DTH  (I,  2)  =  ANGR  (I)  *  0.00125 
Where  I  =  1,  3 

ANGR  =  Angular  Rate 

The  total  angular  increment  fron  the  previous  update  is; 

TDTH  =  DTH  (1,2)  +  DTH  (1,1) 

Where 

DTH  (1,1)  =  Previous  angular  increment,  (the  value  is  always 
less  than  0.0109  Rad.) 

If  TKIH  is  greater  than  0.0109,  set 

\ 

DTH  (1 , 2)  =0 . 0078125*  /  TDTH  \ 

\  ItDTHI;’ 

DTH  (1,1)  =  TDTH  -  DTH  (1,2) 

Otherwise,  set 

DTH  (1,2)  =  0.0 
DTH  (1,1)  =  ODTH 


b) 


Update  the  direction  cosine  matrix  elements  D3(l),  D3(2)  and  D3(3) 


by  the  Crowder-Hession  algorithm. 


D3(2)  =  D3(2)  -  (DTH 
D3(3)  =  D3(3)  +  (DIH 
D3(l)  =  D3(l)  +  (DTH 
D3(3)  =  D3(3)  -  (DTH 
P3(l)  =  D3(l)  -  (DTH 
D3(2)  =  D3(2)  +  (DTH 


(3.2)  *  D3(l)) 

(2.2)  *D3(1)) 

(3.2)  *  D3(2)) 

(1.2)  *  D3(2)) 

(2.2)  *  D3(3)) 

(1.2)  *D3(3)) 
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Ihe  (Xttinands  PC  and  QC  to  move  the  rocket  to  a  new  position 
decided  based  upon  the  value  of  D3  (1) ,  D3  (2)  and  D3  (3) . 


X13  is  calculated  as 

X40  =(31.45  *  THEr)-(31.45  *  X40) 

X41  =  X40  -(90992  *  X13)-(48,26  *  X41) 

X13  =  X41 

The  pitch  angle  X13  is  obtained  through  integration  over  X13. 

e)  The  next  roll  angle  of  the  rocket  is  calculated  as; 

TPHI  =(C76  *  PC)-(C77  *  MGaR(l)) 

• 

X16  =  TPHI  -(100  *  X16) 

TPHI  =  5  *  CX16  +(20  *  X16) ) 

Restrict  angle  TPHI  to  within  +  C5, 

TPHI  =  MAX  (TPHI,  -C6) 

TPHI  =  MIN  (TPHI,  C6) 

• 

The  roll  angle  X14  is  dDtained  through  integration  over  X14  which 
is  calculated  as; 

X35  =(31.45  *  TPHI) -(31.45  *  X35.) 

X36  =  X35  -(90992  *  Xl4)-(48,26  *  X36.) 

X14  =  X36. 

3.6  Subroutine  SEAIMAN 

The  call  to  subroutine  RXKFM  in  subroutine  SEATMAN  has  been  changed 
fron  CALL  RDCKEM  (QT)  to  CALL  RCXIKFM  (QT,  QY,  QDYDT)  . 
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4.0 


Surmiaiy 


CSC  cxxrpleted  ths  modifications  to  the  ICARUS  program  described  in 
Section  3  and  began  preliminary  besting. 


-21- 


